Effective Hamiltonian with holomorphic variables 
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<D • The pure-quantum self-consistent harmonic approximation (PQSCHA) permits to study a quan- 

tum system by means of an effective classical Hamiltonian. In this work the PQSCHA is reformulated 
in terms of the holomorphic variables connected to a set of bosonic operators. The holomorphic 
formulation, based on the olomorphic path integral for the Weyl symbol of the density matrix, 
makes it possible to directly approach general Hamiltonians given in terms of bosonic creation and 
annihilation operators. 

The concept of effective potential in quantum statistical mechanics was introduced by Feynman jjj by means of a 
variational method for the imaginary time path integral. The method was improved by Giachetti and Tognetti Q and 
Feynman and Kleinert 0] using the Lagrangian path integral. Several applications to condensed matter have demon- 
strated its usefulness. H The generalization to Hamiltonian path integrals was performed for treating nonstandard 
Hamiltonians, where there is no separation between a quadratic kinetic energy and a configurational potential, and 
• for which the Feynman-Jensen inequality does not hold in general. An effective Hamiltonian for systems with non- 
l — l \ standard Hamiltonian, as for ex. spin systems, was obtained by means of the pure-quantum self- consistent harmonic 
approximation, (PQSCHA). Q Nevertheless, this generalization can not be applied in a natural way to field-theory 
models which are described in terms of creation and destruction operators: the natural classical-like counterpart of 
. this models requires, in the Bose case, holomorphic variables. We present here the derivation of the effective Hamil- 
tonian in the framework of PQSCHA for Bose systems by using an holomorphic path integral for the Weyl symbol of 
<^> the density matrix p = exp{— f3H}. 

Consider a system of N bosons and let d = {a M } M= i...jv, o, = {aJ,j>=i...A r be the creation and annihilation 
operators. The quantum statistical average of an operator d(a\a) has a natural expression in terms of the Weyl 
0^ , symbols of O and p, respectively 0(z* , z) and p(z* , z), 

g I (6) = \Tr{6p) = | J dp(z*, z) 0(z*, z) p(z\ z) , (1) 

. where Z = exp(—f3F) = Tr(p), the complex variables (z*, z M ) clearly are the Weyl symbols of (d^,a M ), and the 
£^ ' measure is defined as dfi(z*,z) = f\ (dz*dz IJ ,/2Tri). 

The idea is to decompose the path-integral expression for p(z*, z) into a first sum over all paths with the same 
average point defined as ^ Jq du (z* (u), z(u)), and a second sum over average points. For this we introduce in the 
path-integral a resolution of the identity that fixes the average point to (z*, z) and we split the integration over the 
latter, defining the reduced density 
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p(z*,z;z*,z)= J V[z*(u),z(u)]s((z*,z) - i ^ du(z*{u),z(u))^e S ^ z *^' z ^ , 



(2) 



where the Euclidean action S[z*(u), z(u)] is given by 
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S[z*(u),z(u)] = I du 



\ [ t z*{u)z{u)-h*{u)z{u)]-'H{z*{u),z{u)) 



-\ [V(O)*03) - t z*(f3)z(0)] + - [%z*((3) z*(0)]z - ***[*G9) - *(0)]] . (3) 

We take p(z*, z; z*, z) as an unnormalized probability distribution in the variables (z* z) and define its normalization 
constant as exp(~ (3H c e(z*, z)) , so that 
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p(z*, z; 2* z) = e - Wfi («*, *) p(z* » ; 2*, 2) , 

and the average of can then be written 
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dz*dz^ 
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dp,(z*,z)0(z*,z)V(z*,z;z*,z) 



■PH eS (z*, z) 



(4) 
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It is natural to interpret exp(—(3H.cs(z*,z)) as a classical-like effective density, whereas the probability distribution 
V(z*,z;z*,z) describes the particle fluctuations around the point (z* z). In the classical limit it can be seen that 
V(z*, z; z*, z) — > 5 ((z*, z) — (z* z)) and exp(— f3H e g(z*, z)) tends to the classical Boltzmann factor; it follows that the 
probability V describes the pure-quantum fluctuations of the particle thus providing a separation between classical-like 
and pure-quantum contribution to (O) . 

The evaluation of the reduced density p(z*, z; z*, z) will be done in a self-consistent approximation replacing 
H.(z*(u), z(u)) in the action (^) with a trial Hamiltonian quadratic in the displacements from the average point, 
namely 

H (z\ z; z*, z) = %*E(z*, z)£ + \[%F{z*, z)£ + c.c] + w(z*, z) . (6) 

where (£*, £) = (z* — z* z— z) and E = {E^}, F — {F^}, and w are parameters depending on (z* z) that are to be 
optimized. The evaluation of p(z*, z; z* z) is done by diagonalisation of the quadratic term, 
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(7) 



with a canonical transformation (£*, £) that, thanks to the use of path integral with Weyl symbols, preserves 

at the same time the functional measure and the form of the action. 
The explicit result for Vq and Ti. c g is: 



r (z*,z;z*,z)=]l(^^e- 2 ^k/C(f k )^ ^ H cS (z*,z) = w(z*,z) + ±£> 
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fk 



(8) 



where fu{z*,z) = (3wk{z*,z)/2 and £(/&) = coth/fc — f k 1 is the Langevin function. We will denote the TVaverages 
by double brackets ((■}}; of course Pq turns out to be a Gaussian centered in (z* z), so it is fully determined through 
its second moments ((ZtZk')) = Skk 1 £(/fc)/2, = ((£fc£fe')) = 0- The linear canonical transformation made can 

be used to express in a simple way the moments for the original variables, (£«>£/z)- 

In order to determine the parameters of TLq in a self-consistent way we require that the original and the trial 
Hamiltonian, and their second derivatives, have the same "Po-averages: 



({H(z*+c,z+0)) = ({n (z*+e,zH;z*,z))) =J2^k(z*,z) 



+ w(z*, z) . 



(9) 



((d z ,d z ,Mz*+e, = ((d z .d z „Ho(z*+e, zH; **, *))) = E^z*, z) , 

((d^d z ,Mz*+€, z+Z))) = ((dt&Mz'+e, z+£; z* z))) = F^(z*, z) , 

((d z .d zZ H(z*H* = ((d z ,d z ,Ho(z*+e, z+t; z\ z))) = F^z*, z) . (10) 

where ((•}) stands for the Gaussian average given by Vq. 

Finally we observe that a solution of the general problem for arbitrary values of (z* z) is rather difficult to determine 
for many degrees of freedom: a further simplification is in order. This is the low-coupling approximation Q (LCA) 
and its main purpose is to make the averages (( • )) independent of the phase-space point (z* z), so that the above 
self-consistent equations are to be solved only once. The simplest way, consists in expanding the matrices E(z*,z), 
F(z*, z) around a self-consistent minimum of 7i ff • The effective Hamiltonian and the expressions found for thermal 
averages are parallel to those shown in Ref. Q . 

In perspective, one can think of translating this formalism to the fermionic case, getting "classical" expressions in 
terms of Grassman variables. 
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